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^ , Abstract 

' This paper considers the asymptotic behavior of two-source block-fading single-antenna Gaussian 

^ I I interference channels in the high-SNR regime by means of the diversity-multiplexing tradeoff. We 

. . consider a general setting where the users and the average channel gains are not restricted to be symmetric. 

C/3 , 

Our results are not just extensions of previous results for symmetric networks, as our setting covers 
^ ' scenarios that are not possible under the symmetric assumption, such as the case of "mixed" interference, 

' i.e., when difference sources have different distances from their intended receivers. We derive upper and 

00 : 

l/^ , lower bounds on the diversity. We show that for a fairly large set of channel parameters the two bounds 

■ coincides. 

od ■ 

o : 

^ ■ I. Introduction 

O ■ 

K*" ■ Wireless networks deal with two fundamental limits that make the communication problem challenging 

^ ' and interesting. On the one hand, simultaneous communications from uncoordinated users create undesired 

. - - i interference. On the other hand, fluctuations of the channel condition due to multi-path and mobility 

cause signals to fade randomly. In today's cellular and ad-hoc networks orthogonalization techniques, 
such as F/T/C/SDMA, are employed to avoid interference. However, although leading to simple network 
architectures, interference avoidance techniques are suboptimal in terms of achievable rates. Moreover, 
the relative strength of the intended data signal and the interference signals changes over time due 
to fading. This makes fixed channel access strategies suboptimal. Thus, understanding how to deal 
simultaneously with interference and with fading holds the key to the deployment of future broadband 
wireless networks. The simplest model for analyzing these problems jointly is the two-source Block- 
Fading Gaussian InterFerence Channel (BF-GIFC). 
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It is well know that the Han-Kobayashi (HK) [1] scheme with superposition coding, rate splitting, and 
joint decoding, gives the largest known achievable rate region for GIFC without fading. Several outer 
bounds are known in the literature for GIFC without fading [2]-[7]. In particular, Etkin et al. [4] showed 
that a simple rate splitting strategy in the HK scheme is within one bit/sec/Hz of the capacity region of 
Gaussian unfaded GIFCs for any possible channel parameters. In [4], all interfering signals above the 
noise floor are decoded, that is, the private messages -which are treated as noise- are assigned a transmit 
power such that they are going to be received at, or below, the level of the noise. In doing so, roughly 
speaking, the effective noise power at the receiver is at most doubled, thus giving a rate penalty of at 
most 1 bit/sec/Hz. 

Recently, GIFCs with fading were considered in [?], [9]-[15]. 

For ergodic channels, such as fast fading channels, the (Shannon) capacity is the performance measure 
of the ultimate system performance. In [16], it was showed that the sum-rate ergodic capacity of a 
if-source fading GIFC scales Unearly with the number of sources. In [9], the sum-rate capacity of a 
two-source strong ergodic fading GIFCs was shown to be equal to that of the corresponding compound 
MAC. In [10], optimal power allocation poUcies for outer and iimer bounds for ergodic fading GIFCs 
with perfect transmitter CSI were derived. 

For slow fading channels, the proper measurement of performance is the outage capacity. In particular, 
the Diversity Multiplex Tradeoff (DMT) [17], quantifies the tradeoff between rate and outage probability 
as the Signal to Noise Ratio (SNR) grows to infinity. In [11] the DMT of symmetric two-source BF- 
GIFCs was studied based on the "within one bit" outer bound of [4]. The authors of [11] claimed that the 
derived DMT is actually achievable because the "one bit penalty" for using a simple HK strategy vanishes 
at high SNR. However, the achievability of the "within one bit" outer bound requires a very specific rate 
sphtting in the HK achievable scheme that depends on the instantaneous fading values. Hence, as pointed 
out in [I2]-[I4], [18], [19] the DMT derived in [II] is achievable only if the transmitters know the 
instantaneous fading values perfectly. In the the case of no channel state information at the transmitter 
(TXCSI) the DMT of [II] is an upper-bound on the actual DMT. 

The DMT of BF-GIFCs without TXCSI is the subject of investigation of [12]-[15], [18], [19] as well 
as of this work. In [12], it was proved that in strong interference joint decoding of all message at all 
destinations achieves the DMT outer-bound of [11]. In [13], it was showed that multilevel superposition 
coding achieves the DMT of any two-source BF-GIFC; however, no expUcit formula was given for more 
than two levels of superposition. In [15], it was showed that the DMT of BF-GIFCs reduces to that of 
Multiple Access Channel (MAC) if transmitters are not aware of the channel gains. In [14], it was shown 
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that one bit of TXCSI suffices to aciiieve the optimum DMT for certain rages of channel parameters. 

The works [11], [12], [14], [15] focused on two-source symmetric networks, that is to say, networks 
for which the average SNR and the average Interference to Noise Ratio (INR) at all receivers are the 
same. In this work, we consider two-source asymmetric GIFCs as in our conference papers [13], [18], 
[19]. In [18], [19] we generaUzed the DMT outer-bound of [11] to asymmetric networks and studied HK 
achievable schemes with and without rate splitting. It should be point out that our results are not just a 
generalization of the symmetric network results. Our setting covers all possible classes of channels and 
includes channels not possible under the symmetric assumptions, such as the case of "mixed" interference. 
Mixed interference occurs in practice when sources have different distances from their intended receivers 
and is the most practical scenario for wireless networks. 

We assume that the channel variations over time are sufficiently slow so that the fading coefficients 
may be considered as fixed for the whole codeword duration (i.e., block fading assumption). We assume 
that the receivers know perfectly the channel realization, but the transmitters do not. In this case, if 
the instantaneous fading realization is such that the transmission rates cannot be reliably decoded, the 
system is said to experience outage. In an outage setting without TXCSI, it is not clear that a fixed 
rate splitting strategy can actually achieve the DMT upper bound of [11]. Here we consider both HK 
achievable strategies with and without rate splitting. In the case of rate splitting, we consider the case 
where the average received power of the signals that are treated as noise is set below the noise floor, as 
done in [4] for the unfaded case. We also generalized the outer-bound of [11] to asymmetric networks. 
We show that for a very wide range of channel parameters, the inner and outer bound meet. In particular, 
rate splitting improves the achievable DMT in weak and mixed interference channel. 

The rest of the paper is organized as follows: Section |II] presents the system model and the problem 
formulation; Section Hill and ITVl present DMT upper and lower bounds, respecetively; Section W\ presents 
numerical results; Section W\\ concludes the paper. 

II. Channel Model 
A two-source single-antenna Rayleigh fading GIFC in standard form is defined as: 

Yu = HuiXi + Hu2X2 + Zu^C, (1) 

where the noises Zu ~ A/^(0, 1) and the inputs are subject to the average power constraint E[|X„p] < P^, 
u € {1,2}. We assume the channel to be block-fading and that each codeword spans one fading block, 
i.e., no coding across multiple blocks is allowed. Moreover, we assume arbitrarily large block lengths. The 
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receivers are assumed to perfectly know the fading realization {Hii,Hi2, H21, H22), while the transmitters 
are not. In the rest of the paper we parameterize the received SNR/INRs as 

¥.[\H,uXu\'']=n\Hcu?]Pu = x^-, /3cn€M+, (c,n) e {1,2} X {1,2}, (2) 

for some x > 1, and the transmission rates as 

Ru = \og{l + x'-), r„GM+, uG{l,2}. (3) 

We focus our analysis on the high-SNR regime, that is, in the limit for x +00. Notice that, although 
we impose that the channel gains /3's and the rates r's to be non-negative, the results derived in the 
following can be extended to any /3's and r's by replacing each (5 with = max{0,/3} and each r 
with [r]+ = max{0,r}. 

a ) Capacity outer bound: The capacity region of GIFC is not known in general. Recently, Etkin et 
al. [4] proposed a novel outer bound for the capacity region of unfaded GIFC that is shown to be "within 
one bit" of a simplified version of the Han-Kobahyashy [1] achievable region. More precisely, let 

Hi, = ^n\Hij?]x-^^^ ei^-- (4) 

where {jijjOij) are iid for all and with 9ij uniformly distributed on [0, 27r] and independent of 

7jj. In [17] it is shown that in the limit for x —i- +00, the random variables 'jij are asymptotically iid 
negative exponential with mean log(x). By using the parameterization Q, and Q, for each fading 



realization (711,712,721,722)5 the "within one bit" outer bound of [4] can be written 

"T^ETW = |(7ii>7i2,72i,722) e ■ 

log(l + a;''^ < log(l + a;^"-Ti^) 
log(l + x''^) < log(l + x^''-^'') 

^/3ii-7ii 

log(l + x'^) + log(l + x'^) < log(l + ^^—^—^) 

+ log(l + x^^^-T^^ + x^"^""^") 

log(l + x'^)+ log(l + x'^) < log(l + ^-—^—^J 

+ log(l + + x^i^-^i^) 

log(l + x''^) + log(l + a;''") < 

log(l + xf^'^"^'' + 



^./^ii— 711 



1 + a;/^2i-72i 

/^/522— 722 

+ log(l + x^^i-T^i + — o— — ) 

21og(l + + log(l + x'^)< log(l + Y^^^^) 
+ log(l + x^"-'^" + 

f./322— 722 



+ log(l + x^--^- + — ) 



log(l + + 21og(l + x^^) < log(l + ^-—^—^) 

+ log(l + x^''"'^^^" + x^'^"^=') 

7./311-711 "1 
+ log(l +x^^^-^^^ + — — ^ ) }. 

l+a;/321-721^j 
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A. Capacity inner bound 

The HK achievable region, in a fomi that matches the rate bounds of ([5]), can be found in [20], and 
is given by: 

"^HK.complete = U | ("^I'^s) ^ 1^+ : 

P(Q,WuW2,XuX2) ^ 



Ri < I{Xi 

R2 < I{X2 

Ri + R2< I{X2 
Ri + R2< I{Xi 
Ri + R2< I{Xi 
2Ri + R2< I{Xi 
Ri + 2R2 < I{X2 



yi\w2Q); 

y2\WiQ); 

VFi;y2|Q)+/(^i;n|W'it^2,Q); 

VF2;yi|Q)+/(X2;y2|^i^2,Q); 

W2-^Yi\Wi,Q) + I{X2,Wi;Y2\W2,Q); 

W2; Yi\Q) + I{Xi-Yi\Wi,W2, Q) + I{X2,Wi- Y2\W2,Q) 

Wi- Y2\Q) + /(X2; Y2\Wi,W2, Q) + I{Xi,W2; Yi\WuQ) 



The region T^HK.compictc is difficult to evaluate because it requires an optimization with respect to 
the joint distribution P{Q,Wi,W2,Xi, X2), where Q is a time-sharing random variable and (Wi, VF2) 
has the meaning of common information decoded at both receivers. In order to have a region that can 
be evaluated easily, it is customary to assume jointly Gaussian input P{Wi,W2,Xi,X2\Q) without 
time sharing, that is, the random variable Q is a deterministic constant. We set Wu ~ ^{0, Pu,common) 
independent of T„ ~ M{0, Pu,private) and let Xu = Wu + T„ such that the total power constraint is met 
with equality, i.e., Pu = Pu,private + Pu,common, for u G {1, 2}. We further parameterize the ratio of the 
average private power to the total average power for a given user u as 

1 



a,. 



G[0,1], 6„G 



1 + X^" 

so that the average receive SNR/INR's on channel c € {1,2} are 



(6) 



m\H, 



cu\ ]Ru,private 



auE[\H,unPu 



X 



l + x" 



n\Hcu\^]Pu. 



common 



{I - au)E[\Hcu\^]Pu 



l + x'' 



-, uG{l,2}. 



Moreover, following [4], we set bi = (32i and 62 = /9i2 so that the average interfering private power is 
below the level of the noise, that is, ^[\Hcu\'^]Pu,private = irf^ ^ 1- 
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With these choices, 7^HK,compictc reduces to 
T^HK = \ (711,712,721,722) G K 



^/3ii— 711 



log(l + < log 1 + — ^^i^r^ (7a) 



1 + 



log(l+x--)<log(l+ ^ , ....... ) 

\ + l+X + f21 / 

/ ^/322 - 722 I a;''2i-"'2i 

log(l + x^^) + log(l + x'-^) < log ( 1 + . . 



1 + 



l+2:+f2i 



-711 



2;/322 -T22 
J] 

1 + 



X^ll~^ll I X^12~^12 



1 + 



/ 2;'^ll-7ll + X'3l2-712 

2iog(i + x'-o + iog(i + x'-^) < log 1 + ^ ..i.'-tr '" - 

\ + l + X + fl2 




(7b) 



iog(i + x'-^) + iog(i + x'-^) < log 1 + ^ ,,i.i!r; '" - 

V + l + X + ft2 / 



+ log|l + T^^?5^ (Vd) 



log(l + x'-^) + log(l + x'-^)<log 1+ . 

\ 1- + l+X + '3l2 / 

/ x^22-T22 3;/321-721 \ 

,1 (11 l + X + ft2 + l + x-''21 \ ,„ , 

+ log 1 + . . ^..,-.2, (7e) 



(7f) 



/ 2;'^22-722 I X'321-T21 \ 

log(l + x'-O + 21og(l + x^-^-) < log 1 + ^ 

V l+X + f21 J 

(^/322 -722 \ / 3;/3l2-"^12 \ >^ 

1 + , ..21-.21 + log 1 + - . (7g) 

i + I+X + P2I ) V ^ + l+X + fl2 )> 



B. Diversity 

The probability of outage Pout (''i , ^^2 ) is defined as the probability that the fading realization (711 , 712 , 721 , 722 ) 
is such that the rate pair (ri, r2) cannot be decoded. By using the outer bound region T^etw in © and 



the inner bound region T^hk in © we can bound the outage probabiUty as 

1 - P[(7ll,7l2,721,722) G ^ETW] < Fout(?'l,T2) < 1 - IP[(7ll , 7l2 , 721 , 722) GT^Hk]- 

The diversity, or the high-SNR exponent of the outage probability, is defined as 

,( ^ r -log(Pout(r-i,r2)) 
d(ri,r2)= hrn -— , 

and it is bounded by 

dBK{ri,r2) < d{ri,r2) < dETw(n,'r2), (8) 

where dETw{^i,f2) and (^hkI'^i, ^^2) are defined similarly to d{ri,r2). 

The rest of the paper is devoted to the evaluation of dETw{fi:i"2) and duK{ri,r2). 

III. Diversity upper bound 
By using the Laplace's integration method as in [17] we obtain 

dETw{ri,r2) = mill {711 + 712 + 721 + 722} (9) 
where T^etw is the large-x approximation of T^etw in © and is given by 

^ETW = I (711, 712, 721, 722) G ■■ Xij = [Pij - 7ij] + , 

ri < Xn (10a) 
r2 < X22 (10b) 

rs = ri + r2 < [Xn - X2i]+ + max{X2i, X22} 

= maxjXii, X21} + max{X22, -'^21} - -^21 (10c) 
Ts = ri + r2 < [X22 - Xu]'^ + max{Xi2, Xn} 

= max{Xii,Xi2} + max{X22,Xi2} - X12 (lOd) 
rs = ri + r2 < max{Xi2, Xn - X21} + max{X2i , X22 - -^^12} 

= max{Xii, X21 + X12} + max{X22, X21 + X12} - {X21 + X12) (lOe) 
rf = 2ri + r2 < [Xn - X2i]'^ + max{Xii,Xi2} + max{X2i, X22 - Xu} 

= max{Xii,X2i} + max{Xii,Xi2} + max{X22, ^21 + ^12} - (^21 - ^12) (lOf) 
rg = ri + 2r2 < [X22 - ^12]+ + max{X22, X21} + max{Xi2, Xn - X21} 

= max{X22,X2i} + max{X22,Xi2} + max{Xii , X21 + Xu} - {X21 - X^)^ (lOg) 
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where = max{0, x}. 

The optimization problem in ^ can be solved as follows: since the complement of T^etw is the union 
of the complement of the conditions (llOal ) through ( 10g[ ), by applying the union bound as in [21] it can 
be shown that the diversity in ^ evaluates to 



max {-^11 

X's do NOT satisfy equation (TOt ) 



Xi2 + X21 + X22}- 



d^Q^) = Pn + P12 + P21 + P22 
We have: 

• The diversity ti |yQ^ (corresponding to the constraint (llOal l) is: 

d []M =/5ii - max{Xii} 

subj. to < Xu < Pii, Xii < ri, 
=/3ii - min{^ii,ri} = -min{0,ri - Pu} = max{0,/3ii - n} 
= [/3ii-ri]+. 

• Similarly to (i |jQ^ , the diversity d|[Job} (corresponding to the constraint (llObl l) is: 

^ ITObt = [/^22 - r2]+. 

• The diversity d|JOcl (corresponding to the constraint (llOcI l) is: 

d^Q^ =Pii + P21 + P22 - max{Xii + X21 + X22} 

subj. to < Xn < Pu, < X21 < P21, 0<X22< P22, 
and to maxjXii, X21} + max{X22 , -^^21 } - X21 < Vg = ri + r2. 
We start by re-writing the last constraint as follows: 

X + Z<rs+Y 



max{X, Y} + max{Z, ¥} <rs + Y 



Y + Z<r, + Y 
X + Y <rs + Y 

Y + Y <rs + Y 



X + Z<rs+Y 
max{X, y, Z} < rg 
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which implies 

max{X2i + (Xii +X22)} 

subj. to < X21 < min{rs,/32i}, < Xn < mm{rs,Pii}, < X22 < mm{rs, f322}, 

and to Xii + X22 < + X21, 
= max|x2i + minj min{rs,/3ii} + mm{rs, P22},rs + ^21^^ 

subj. to < X21 < min{rs,/32i}, 
= min| min{rs,/5ii} + min{rs,/322} + min{rs, /52i}, + 2min{rs,/32i}|. 
Hence we obtain: 

- minj min{rs,/3ii} + minjrs, /322} + minjrs, /32i}, + 2 min{rs, /32i}| 

= max{[/3ii - r,]+ + [(322 - rs]+ + [P21 - rs]+,Pn + P22 - 2r, + |/32i - r,|}. 

Remark: In the symmetric case, with f3ii = (322 = 1> P12 = /?2i = a and ri = r2 = r (i.e., Ts = 2r), 
d|YQ^ reduces to 

dm,sym = 2 max { [1 - 2r]+ + [| - r] + , (1 - 2r) + 

= 2 ( [A]+ + [B]^ + \[B]- - [AY 



a 

r 

2 



where 



\A=l-2r,B= 



X = [x\^ — [x] : = max{0, x] > 0, [x] = — min{0, x} > 0. Vx G M. 



The corresponding bound in [11] is 



4 =2[1 - r - min(r, ^)]"'" + [a - 2r] + 

= 2 f [^]+ + [5]+ + \[A - B]+ - [A]+ + [B]- 



\A=l-2r,B=^-r, 



which can be easily shown to be equivalent to d jjQ^ since for ^ > : — 5]+ — [^]+ + 
[B]+ = A- B - [A]+ + [B]+ = -[A]- + [B]-, and A < B : [A- B]+ - [A]+ + [B]+ = 
-[A]+ + [B]+ <0. 

The diversity c?|Xod| (corresponding to the constraint (llOdl l) is as (i |jQ^ but with /32i replaced by 
j3i2, i.e., with the role of the users swapped. 
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Remark: In the symmetric case, d^Q^^^^ = d^^^^. 

The diversity (i jjQg| (corresponding to the constraint (llOel l) is as ti |jQ^ but with P21 + /?i2 instead 
of P21. 

Remark: In the symmetric case, c? |Yq^ coincides with dc in [11]. 
The diversity c^Uofl (corresponding to the constraint (llOfl l) is: 



rf |TOft =/3ii + P12 + 021 + (322 - max{X + Y + Z + W} 

subj. to < X < 0<W < /3i2, < y < /321, < Z < /522, 

and to max{X, y} + max{X, VF} + max{Z, Y + W} - {Y + W) < rj = 2ri + r2. 

The last constraint can be rewritten as 

{X -Y) + {X -W) <rf - max{Z, Y + W} 
X -Y <rf- max{Z, Y + W} 
X -W <rf- max{Z, y + T^} 
< r/ - max{Z,y + W} 
that together with X < /Jn gives 



2X < rj - [Z - (y + VF)]- 
max{Z,y + W} <rf 



maxjA ) = mm < pn, — 

Next, the function 

max{X} + {Y + W) = min "^-^ ~ ~ ^ } + (X + W), 



is increasing in (Y + W), and since the constraints 

W < f3u, Y < P21, Z < P22, max{Z,y + W} < rj, 

are equivalent to 

Y + W < min{/52i +/3i2,r/}, Z < min{/?22, r/}, 

we have 

max{X + y + 1^} = mm {/3n /^ " " + ^"^^^ } + min{r„ /32i + /3i2}. 

Finally, the function 

Z + max{X + y + ly} = mm {/3n, " " + ^"^^^ } + mm{r„ /32i + A2} + Z, 
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is also increasing in Z, hence, subject to Z < mm{rf, P22}, we finally have 

(i fTOft = - max{X + Y + W + Z} + {(3u + /3i2 + P21 + (^22) 

- miTi//^ ^/ - [minjr/, /322} - min{r/, /^ai + /?i2}]+ ] r 0,01 „,.^r^ ^ n 

— — mm < pii, > — mm|rj, P21 + P12 j — min|rj, /:)22| 

+ iPii + /3i2 + P21 + P22) 

Tf -[(322 - {f32i + f3i2)+a-b]- 



+ a + b, a= [ip2i + (3i2) - b = [P22 - 



2 

Remark: In the symmetric case, (j |jQf| reduces to 



,sym 2 

which is not equivalent to in [11]. In fact, it turns out that dd in [11] is not correct. Consider 
the following numerical example: let ri = r2 = 0.4, /3u = P22 = 1> P12 = P21 = 0.5, which 
corresponds to rj = 1.2, a = 0.5 in [11]. The optimization problem for dd is 

dd = min{7ii + 712 + 721 + 722} 
subj.to [[1 - 7ii]+ - [a- 712]^]^ + max([l - 711]+, [a - 721]+) 
+ max([a - 712]^, [1 - 722]^ - [a - 721] + ) < r/ 

It can be easily verified that 711 = 0.4, 712 = 721 = 722 = is a feasible solution that gives 
(711 + 712 + 721, 722) = 0.4. However, according 

dd = max{(l + (1 — 3r)^ + (2a — Sr)"*", min{[3 — 3r — min(3r, 2q)]''', max(l, 2 — 3r — min(3r, 2a))} 

we have dd = 0.8. 

Finally di^Ogl (corresponding to the constraint ( 10g| )) is as d |jQf| but with (^22 instead of Pu and 



Vg instead of rj, i.e., the role of the users is swapped. 

IV. Diversity lower bound 

The evaluation of the diversity lower bound dux in ^ can be carried out similarly to the evaluation 
of the diversity upper bound c^etw in the previous section. We will consider both the case of no rate 
splitting and the particular choice of the power split among common and private messages inspired by [4] 
which led to 
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A. Diversity lower bound without rate splitting 

Without rate splitting in tlie HK region, a user either sends all private information or all common 
information. These two modes of operation correspond to either treating the interference as noise at the 
receiver, or performing joint decoding as in a MAC channel. 

Consider first the case where the interference is treated as noise. User 1 can be successfully decoded 
at receiver 1 by treating user 2 as noise if 

{(711,712) G m2 : log(l + x^^)< log + j^^^) }• 

By following the same approach used in the derivation of the diversity upper bound, we have that the 
exponent of the probability that user 1 cannot be decoded successfully at receiver 1 by treating user 2 

as noise is given by: 

dmi = Pii + (3i2 - max{Xii + X12} 

subj. to < Xu < /3ii, < X12 < P12, [Xu - Xi2]+ < n 
=/3ii + P12 - max{Xii + Xu} 

subj. to < Xn <Pn, 0< ^12 < P12, Xu < X^ + n 
=Pii + P12 - max{min{/3ii, X12 + n} + Xu} 

subj. to < X12 < P12, 
=/3ii + 012 - max{min{/?ii,/?i2 + ri} + ^12} 
= [Pii - ri - /3i2]+, 

Similarly, the exponent of the probabihty that user 2 cannot be successfully decoded at receiver 2 by 
treating user 1 as noise is 

C?NI2 = \P22 -r2- /?2l]"^- 

Consider now the case where the users are jointly decoded. User 1 and user 2 can be successfully 
jointly decoded at receiver 2 as in a MAC if 

{(721,722) G : 

log(l +X''')< log(l + x^^i-T^i) 

log(l + x''^) < log(l + x^''-^'') 

log(l + x''') + log(l + x''^) < log(l + x'^^'-^^' + 0;^=^=-^^=)}. 
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The exponent of the probability that both users cannot be jointly decoded is given by 

dMAC2 = min{[/32i - ri]+, [/?22 - ^2]"^, [/322 - fs]"*" + [1321 - Ts]'^}, rs = n+ r2, 
where the last argument of the minimum in dMAC2 can be derived as follows: 
1^22 + /321 - max{X2i + X22} 

subj. to < X22 < P22, < X21 < P21, max{J^22,^2i} < r^, 
=022 + P21 - min{/32i, Ts} - min{^22, rs} 
=[^22-r,]+ + [;32i-r,]+. 

Similarly, the exponent of the probability that user 1 and user 2 cannot be successfully jointly decoded 
at receiver 1 is 

c^MACi = min{[/3ii - ri]+, [/3i2 - r2]+, [f3u - r^j^ + - rs]~^}, = ri + r2. 
Hence, without rate splitting, we have 

c^HK-wors = max{doo, doi, dio, cili}, (11) 

where "wors" stands for "without rate splitting" and where 

• dii is the diversity when both sources send only private information (which is sum-rate optimal for 
very weak interference unfaded GIFCs [5]) given by 

du = min{dNii,(iNi2}- 

• dio (and similarly for doi but the role of he users swapped) is the diversity when user 1 sends only 
private information and the user 2 sends only common information (which is sum-rate optimal for 
mixed interference unfaded GIFCs [22]) given by 

dio = min{cZNii,dMAC2}- 

• doo is the diversity when both sources send common information (which is optimal for strong 
interference unfaded GIFCs [?]) given by 

doo = min{dMACi,cfMAC2}- 
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Since max{min{a, 61}, min{a, 62}} = min{a, max{6i, 62}}> we further rewrite dnK-wors in (fTTl ) as 

^HK-wors = max { min{(iNii, dNl2}, min{(iNli, dMAC2}, min{dMACi, (^nk}, min{(iMACi, c^mAC2}} 

= max { min{(iNii, max{(iNi2, dMAC2}}, min{dMACi, niax{dNi2, dMAC2}}} 

= min { max{(iNii, dMACi}, max{dNl2, 'iMAC2}} 

= min { dwrp 1 , ^worp 2 } , dworp u = maxj^Ni u , duAC u } , 

that is, the diversity dnK-wors has the following intuitive explanation: each user u G {1,2} chooses the 
best strategy between treating the interference as noise (dmu) and joint decoding (dMACu), which gives 
diversity d^^pu, and the overall diversity is dominated by the the worst user. 

B. Diversity lower bound with rate splitting 

In the general case of rate splitting, with the power split indicated in Section [HI we have: 

dHK(n,?^2)= min {711+712+721+722} (12) 

where TZhk is the large-x approximation of TZuk in © and is given by: 

^HK = I (711 1 712 , 721 , 722) G : Xij = Pij - Jij , 

n < [Xn]+ (13a) 

r2 < [^22]+ (13b) 

rs = n + r2< [max{X22,X2i}]+ + [Xu - /32i]+ (13c) 

rs = n + r2< [max{Xn,Xi2}]+ + [X22 - /3i2]+ (13d) 

rs = n + r2< [max{Xn - /32i, ^12}]+ + [max{X22 - /9i2, ^21}]+ (13e) 
rj ^ 2ri + r2 < [max{Xn, Xi2}]+ + [Xu - /?2i]+ 

+ [max{X22-/3i2,X2i}]+ (13f) 
rg = n + 2r2 < [max{X22, ^21}]+ + [X22 - Pu]^ 

+ [max{Xn-/32i,Xi2}]+|. (13g) 
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We finally have that the diversity is given by: 
dHK = mill {du3j)}, 

z — a. ..g 

drm. = Pn + Pi2 + P21 + P22 - max {Xn + Xu + X21 + X22}. 

^^^> X's do NOT satisfy equation fB]^) 

where 

• The diversity d|Y3^ and ^ ^3^ (corresponding to the constraint (I13al ) and (I13bl ). respectively) are: 
and 

^IT3bt = ^ITObt = [/?22-r2]+, 

as for the upper bound. 

• The diversity d|j33| (corresponding to the constraint (I13cl )) is: 

rf |T33t =/?22 + /321 + /?ii - max{X22 + X21 + Xn} 
subj. to X22 < P22, X21 < P21, Xu < pn, 
and to [max{X22, ^21}]+ + [^ii - ^21]^ < = ri + r2. 
If rewrite the optimization domain as 

-^^22 < P22, X21 < P21, max{X22,X2i} < ars, 
Xn <Pii, Xn-I32i < (1 - a)r„ 
QG [0, 1], 

we immediately obtain 

^njcl = /^22 + P21 + Pu - max \ min{/?22, arg} + min{/?2i, arj + min{/3ii, (1 - a)rs + P21} \ 
' ' oe[o,i] J 

= mill {[P22 - ar,]+ + [/?2i - ar,]+ + [pu - P21 - (1 - a)r«] + }, 

aG[0,l] 

and the optimal a (see Appendix lAl) is 

« |l33t = min{l,max{/322,/32i}As} • 

Remark: In general, it is difficult to compare the lower bound c^U^^J with the upper bound c^|Yoc|- 
By inspection, if > P21 > max{/3ii, /322} the two bounds meet. It is possible that this condition 
we can be relaxed. 
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The diversity d|Y3d| (corresponding to the constraint (I13db ) is as d iJg^ but with the role of the 
users reversed. 

The diversity d jjg^ (corresponding to the constraint ( |13e| )) is: 

d |T33t =P22 + /321 + + /3i2 - max{X22 + X21 + Xn + ^12} 
subj. to X22 < /522, X21 < (321, Xn < Pn, ^12 < /5i2, 
and to [max{Xii - /32i,Xi2}]~^ + [max{X22 - /3i2, -^^21}]"^ < r^- 
If rewrite the optimization domain as 

-^^22 < 1^22, X21 < P21, max{X22 - ^12,-^^21} < ars, 

Xu < Pii, X12 < P12, max{Xii - /?2i,^i2} < (1 - a)rs, 

a G [0,1], 

we immediately obtain 

rf |T33t = P22 + P21 + Pn + P12 

— max <^ min{/322, /3i2 + ars} + min{/?2i, ar^} + min{/?ii, /32i + (1 - o:)rs} + min{/3i2, (1 

og[o,i] 

= min \[p22 - P12 - ar,]+ + [P21 - ar,]+ + [Pu - P21 - (1 - a)rs]+ + [/3i2 - (1 - a)r,] 

and the optimal a (see Appendix iBt is 

A . f max{^22 - /32i,/?2i} 1 
«{l3i = ™in|l' 1- 

The diversity (corresponding to the constraint (I13fl )) is: 

d |T3ft =/?22 + P21 + + P12 - max{X22 + X21 + Xu + Xu} 
subj. to X22 < P22, X21 < P21, Xu < Pu, X12 < P12, 
and to [max{Xn, Xi2}]+ + [Xn - p2i]^ + [max{X22 - /3i2, ^21}]+ < r/ 
If rewrite the optimization domain as 

^22 < P22, X21 < P21, max{X22 - /3i2,X2i} < arj, 

Xu < Pu, X12 < P12, max{Xii,Xi2} + [Xu - P21V < (1 " a)rf, 

oe [0,1], 
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and max{Xii,Xi2} + [Xu - /32i]+ < (1 - o.)rf as 

^11 + ^11 < P21 + (1 - a)r/ 
^ii+/?2i </?2i + (l-a)r/ 
^12 + ^11 < P21 + (1 - a)r/ 
^12 + /321 < P21 + (1 - a)r/ 
we obtain that the optimization domain is 

X22 < min{/322,/5i2 + arj}, X21 < min{/?2i, ar/}, 

P21 ~l" (1 — f 

Xn < min{/3ii, (1 - a)r/, Xu < min{/?i2, (1 - a)r/}, 

^11 + ^12 < P21 + (1 - a)r/, a e [0, 1], 
and we immediately obtain 

rf| [T3ft = ^22 + /321 + /3ii + P12 - max | min{/322 , /?12 + ar/} + min{/32i, ar/} 

+ min|/32i + (1 - a)rf, mm{Pn, (1 - a)?^/, — ^^iL/j _^ min{/3i2, (1 - a)r/}||, 

and the optimal a can be found in Appendix O 
• The diversity dnj^ (corresponding to the constraint (13gl) is as rf jjj^ but with the role of the 
users reversed. 

Remark: In this section we derived the achievable diversity for the case of power split such that the 
average power of the private message at the non-intended receiver is below the noise floor. In Appendix ICl 
we derive the achievable diversity for a generic power split. The expression is quite complex and not very 
insightful. By numerical optimization we found that the particular power split we chose in Section |II] is 
optimal, or very close to optimal, for a very large set of channel parameters {(3u, P12, P21, 1322)- 

V. Numerical results 

In this section we present numerical evaluations of the diversity upper bound c^etw in ® and the 
diversity lower bound without rate splitting dnK-wors in (fTTI ) and the diversity lower bound with rate 
splitting d}iK in (fT2l ) for different values of the channel parameters /9i2, /32i5 /322)- 
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Multiplexing Gain: ^^=^2^' 



Fig. 1. Symmetric channel in weak interference: f3ii — P22 = l,/3i2 = /32i = 0.2. 



A. Symmetric channels 

We first consider symmetric channels. We set the average received power of the direct links to = 
/322 = 1 and the average received power of the cross links to /3i2 = /32i = /? > 0. In Figs. [TJ [2j [3 IH [51 
and[6]we plot the diversity vs. the common multiplexing gain ri = r2 = r for /? = 0.2, 0.5, 0.6, 2/3, 1.1, 
and 1.5, respectively. 

In Figs. [T] and |2] (weak interference) c^etw = c^hk and the dominant constraint at low rate is the 
"single user diversity" d|JOa| ~ ^ ITObt ' while at high rate is the "sum-rate diversity" d |yQg| . 

In Figs. [3] and |4] (weak interference) c^etw = c^hk and the dominant constraint at low rate is the 
"single user diversity"d|jQ^ = d^^, at medium rate is = dr^Q^ and d^^ = d^^, while at 

high rate is c?|Xoa| ~ '^{TObJ again. These figures show that the expression given for dsTW in [11] given 
by c^AL is not correct, as pointed out in a remark earlier. 

In Fig. [5] (strong interference) c^etw 7^ c^hk and ^HK-wors = ^hk- The dominant constraint at low 
rate is the "single user diversity" d|Yoa| ~ ^ fTObt ' medium rate is the "sum-rate diversity" d|Yoc| 
while at high rate is again "single user diversity" d|Xoa| ~ '^(TOb}- 

In Fig. [6l (strong interference) c^etw = c?HK-wors = ^hk and the dominant constraint at low rate is the 
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Qi ' ' ' ' ^- I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Multiplexing Gain: r.|=r2='' 



Fig. 2. Symmetric channel in weak interference: diversity vs. ri = r2 = r for /3ii = P22 ~ l,/3i2 = /32i ~ 0.5. 

"single user diversity" (i |jQ^ = ii|Yob|' while at high rate is the "sum-rate diversity" (i jjQg| . In this case 
no-rate splitting is optimal. In [12] it was show rate splitting is not needed in very strong interference, 
that is, for (3i2 = P21 = P > 2. Here we show numerically that the threshold of 2 for the average 
interference power can be lowered. We found by simulation that c?etw is not achievable for symmetric 
channels for /3 G [0.680, 1.500]. 

B. Asymmetric channels 

In Figs.|7J[8]and|9]we plot the diversity vs. the common multiplexing gain ri = r2 = r for asymmetric 
channels with Pu = /322 = 1 and (3i2 7^ (32i- 

In Fig. |7] (weak interference) c^etw = duK and the dominant constraint at low rate is the "single user 
diversity" d jjQ^ = (j |YQ^ since Pn = ^22, at medium rate is and (j |jQg| , while at high rate is 

d jim = ^{lOb} again. 

In Fig. [8] (mixed interference) c^etw = duK and the dominant constraint at low rate is the "single user 
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Fig. 3. Symmetric channel in weak interference: diversity vs. ri = r2 = r for /3ii = /322 ~ l,/3i2 = P21 ~ 0.6. 



diversity" d^Q^ = d^^, at medium rate is d^Q^ and d^Q^, while at high rate is d^Q^ = d ^ 
again. 

In Fig. |9] (very strong interference) c^etw = ^HK-wors and the dominant constraint is (j |jQ^ = (j |YQ^ . 
In very strong interference, i.e., min{/3i2, /32i} > + (322, the interference is so strong that each user 
can completely remove the unintended signal before decoding its own signal. In this case the capacity 
region, and hence the diversity, is the cartesian product of the single user capacities without interference. 

VI. Conclusion 

In this paper, we analyzed the diversity-multiplexing trade-off of two-source block-fading Gaussian 
interference channels without channel state information at the transmitters. As opposed to previous works, 
we considered generic asymmetric networks. We found that, a simple inner bound based on the HK 
scheme with fixed power split achieves the outer bound based on perfect channel state information at the 
transmitter for wide range of channel parameters. 
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Fig. 5. Symmetric channel in strong interference: diversity vs. ri — r2 = r for /3ii = (322 = 1, 13i2 ~ P21 ~ l-l- 




Multiplexing Gain: '^='2=' 



Fig. 6. Symmetric channel in strong interference: diversity vs. ri — r2 = r for /3ii = (322 = 1, /3i2 ~ P21 ~ 1-5. 
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Multiplexing Gain: '■^=r2='' 



Fig. 7. Asymmetric channel in weak interference: diversity vs. n = r2 = r for /3ii — P22 = l,/3i2 ~ 0.9, /32i = 0.2. 
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Fig. 8. Asymmetric channel with mixed interference: diversity vs. ri = r2 — r for /3ii — P22 ~ l,/3i2 = 1.2, /32i ~ 0.5. 
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Fig. 9. Asymmetric channel in strong interference: diversity vs. n = r2 — r for Pn = (322 ~ l,/3i2 = 3,/32i — 5. 
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Appendix 

A. Optimization of d gj^ 

The optimization problem with respect to a in d jjg^ involves a function of the type 

(i|Y33j(a) = [m- q]+ + [M - a]+ + [-m2 + a]+, m < M. 

For Q < m 

d |X33t (a) = m + M - 2a + [-m2 + a]"^ 

which is decreasing in q since, depending on the value of m2, d. ^ij^ {a) is a straight line of slope -2 or 
-1. For m < a < M 

rf |X33| i(a) = M - a + [-m2 + a]+ 

which is non-increasing in a since, depending on the value of m2, fdct) is a straight line of slope -1 or 
0. For a > M 

^n3ct(«) = [-"i2 + a]+ 

which is non-decreasing in a since, depending on the value of m2, d ^Y^ {a) is a straight line of slope 
or -1-1. Hence the function has a minimum (maybe not unique) at a = M. 

In general, it can be easily shown that the function d ijg^ (a) is fiat in the interval with extreme points 
max{m,m2} and M. This means that the minimum of d ^i^ {a) over a € [0, 1] is achieved by any a 
in the interval with extreme points min{l, max{m, 771,2}} and minjl, Af}. In particular, we can chose 
a = minjl, M}. 

B. Optimization of d^j^ 

The optimization problem with respect to a in rf |Y3^ involves a function of the type 

with TTT-i < Ml and with m2 < M2. With similar reasoning as in Appendix |Al it can be shown that 
the function is flat in the interval with extreme points max{mi, 777,2} and min{Mi,M2}. This means 
that the minimum of d^^{a) over a G [0, 1] is achieved by any a in the interval with extreme points 
min{l, max{777i, 7712}} and min{l, Mi, M2}. In particular, we can chose a = min{l, Mi, M2}. Notice 
that d^2^ = a!|X33j for M2 = +00. 
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C. General rate splitting in the HK achievable region 

The large-x approximation of the HK region for a general power split 
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IE[|-ffcu| ]Pu,private — Cku IE [ 1 | ]Pu 
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l + x'' 
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(711, 712, 721, 722) G ■■ Xij = Pij - jij, 

ri < [[Xn]+ - [X12 - [-fc2]+]+]+ 
r2 < [[^22]+ - [^21 - [-foi]+]+] + 
n + ra < [max{[X22] + , [X21 - [+61]+] + } - [X21 



-61]- 
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n + r2 < [max{[Xn]+, [X12 - [+62]+]+} - [X 
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n + r2 < [max{Xn - [-6i]+,Xi2 - [+62]+} - [X^ - [-62]+]^ 
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^11 - [+61]+ - [^12 - [+62]+]+]" 
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:X22 - [+62]+ - [^21 - [+61]+]+]" 



max{X22 - [-62] + ,X2i - [+61] + } - [X21 
+ r2 < [max{Xii,Xi2 - [+62] + } - [^12 - 
Xu - [-61]+ ~ [X12 - [-62]+]+] + 
max{X22 - [-62] + ,X2i - [+61] + } - [X21 



ri + 2r2 < [max{X22,X2i - [+61] + } - [X21 



-b2V 
-h] 

-hV 



X22 - [-62]+ - [X21 - [-61]+]+]+ 
max{Xn - [-bi]+,Xu - [+62]+} - [X^ 



(14a) 
(14b) 

(14c) 

(14d) 

(14e) 



(14f) 



(14g) 



The evaluation of the diversity can be carried out similarly to the evaluation of the diversity upper and 
lower bounds as done previously. In particular: 

For il4ai (and similarly for (I14bl l): is as "treating the interference as noise" but with interference 
level with [P12 — [+^2]^]^ instead of P12, that is, 

= [^11 - - [+^2]+]+ - ri]+ G [[Pu - P12 - + , [Pu - n]+] 
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The minimum value of d|j4^ is attained for 62 < 0, while the maximum value is attained for &2 > /3i2- 
Recall, 62 = P12 is the power split we chose in the main section of this paper. 

Before the derivation of f^UJ^} = dc (and similarly for rf |Y4dl |)» ^([l4e} ^ ^^'^ '^{I4fl ^ ^^^'^ 



similarly for (in"4^) corresponding to the constraint in (I14cl ). (I14el ). (I14fl) let 

A = [/3n]+, i? = - [61]+] + , 
C = [/?22]+, ^) = [/322 - [62]+]+, 

= [/3i2 - [-62]+]+, F = [/?12 - [62]+] + , 
G = [P21 - [-h]+] + , H = [1321 - [61]+] + , 
Ts = [r-i]+ + [r2] + , rj = 2[ri]+ + [r2] + , = [n] + 2[r2] + . 
For (I14cl l (and similarly for (I14dl l) we need to solve: 
dc = min{7ii + 721 + 712 + 722} 
subj. tomax{C - 722 - [H - 721]+, [G - 721]+ - [H - 721]+}+ + [B- 711 - [F - 712]+]+ < r,. 
We divide the optimization into two steps. First we solve 

del = min{7ii + 721 + 712 + 722} 

subj. tomax{C7 - 722 - [H - 721]+, [G - 721]+ - [H - 721]+}+ < 

and 

dc2 = min{7ii + 721 + 712 + 722} 

subj. to[B - 711 -[F - 712]^]^ < rs2- 

Then we solve 



dc = min{dci + 42 } 

subj. to Tsi + rs2 = rs = ri+ r2. 

The optimization problem has the following three forms: 
1) CASE I: G > H and r, > G - [H]+. 
a) If r,i<G- [H]+, then: 

dc=[C- r,i]+ + [G- r,i]+ + [B - - r.a] 
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b) If min{G -[H]+,rs}< mm{C, G}, then: 

mindc = min{(C - r,i) + (G - r,i) + [B - [F]+ - rs2]+}. 

Increasing one unit of Vgi results in two units decrease of the object function, while increase 
one unit of results in one unit (or less) decrease of the object function. Thus r^i = 
min{G — [i?]"*", Tg}, and rs2 = and 

de=[C-G]+ + [B-[F]+. 

c) If min{G -[H]+,rs}> min{C, G}, then: 

del = mm{[C - r,i]+ + [G - r,i]+ + [B - [F]+ - r,2]+} 

Reasoning as before, rgi should be as large as min{G, G}, thus let's assume Vgi = min{G, G}+ 
r'gi, then the object function turns out to be: 

del = min{[max{G, G} - min{G, G} - r',^]-^ + [B - [F]+ - r^2]+} 

trivially solved as: 

del = [max{G, G} - min{G, G} + [B - [F]+]+ - r;]+. 
Hence, by defining Ti = min{G, G,G — rs}+, and T2 = rg — Ti we obtain 

dei=[C + G- 2Ti + [B- F]+ - T2]+. 

d) If rsi>G- H, then 

de2 = min[G -H- r,i]+ + [B-F- r,2]+ 
= [[C-H] + [B-F]+-rsV 

and hence 

de = min{dci,c?c2}- 

2) CASE 2: G > and < G - 

If Tsi < G — H, by applying the reasoning in the previous subsection and by define Ti = 

min{G, G,G — H, r^}^, and T2 =rs — Ti we obtain 

de=[C + G-2Ti + [B- F]+ - T2]^. 

3) CASE 3:G<H. 



In this case 

dc = [C-H- r,i]+ + [B-F- r,2]+ 

and we trivially have 

m.mdc=[C - H +[B - F - r,]. 

For dl4el l we need to solve: 

de = min{7ii + 721 + 712 + 722} 

subj. to max{i? - 711 - [F - 712]"^, -E - 712 - [F - 712]^}^ 

+ max{D - 722 -[H - 721]^, G - 721 - [-ff - 721]^}^ < 1 
We divide the optimization into two steps. First we solve 

del = min{7ii + 721 + 712 + 722} 

subj. to max{Z) - 722 - [H - 721]"*", G - 721 - [i? - 721]'''}''' < f 

and 

de2 = min{7ii + 721 + 712 + 722} 

subj. tomaxji? - 711 - [F - 712]"'', E - ^12 - [F - 712]^}^ < r 

Then we solve 

de = minjdei + ^2} 
subj. to rsi + rs2 = rs- 

The optimization problem has the following four forms: 
1) CASE 1: G > H and E > F. 

a) If Vsi > G — H and rs2 > E — F, which requires 

rs>[G-H]+ + [E-F]+, 

then 

del = [[D-H]+ + [B-F]+ -rs]. 

b) If Vsi > G - H and rs2 < E - F, which require > [G - H]'^, then 

de2 = [D-H- r,i]+ + [B- r,2]+ + [E - r,2] + 
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and 

minde2 = [B + E-2Ti + [D- H]+ - T2]+ 

where 

Ti = mm{B, E,E-F, rj 
T2 = rs- Ti 



c) If r^i < G — and rs2> E — F, which requires rg > [E — then 

de3 = [£> - + [G - r,i]+ + [B-F- r,2]+ 

and 

de3 = P + G - 2Ti + [B- F]+ - T2] 

where 

Ti = min{L>, G, [G-ii-]+,r,} 
T2 = rs- Ti 

d) If < G - and rs2 < E - F, which requires < [G - + - then 

^4 = [D + G - 2Ti + B + - - r3]+ 

where 

Ti = min{D, G,rs,G-iJ} 

T2 = mm{B, E,rs-Ti,E-F} 
Ts = rs-Ti- T2 

For the four cases above, we have 

de = mm{del,de2,de3, de4}- 

2) CASE 2: G > and < F. 

a) If Tsi > G — H, wiiich requires rg > G — H, then 

dei = [D-H- r,i]+ + [B-F- r,2]+ 
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and 

mindei = [[D - H]+ + [B- F]+ - rs]+ 

b) lfrsi<G-H than 

de2 = [D- r,i]+ + [G- r,i]+ + [B - F - r,2]+ 

and 

de2 = [D + G- 2Ti + [B- F]+ - T2]+ 

where 

Ti = mm{D,G,G - H,rs} 
T2 = rf- Ti 

For these two cases we have 

de = max{dei,de2} 

3) CASE 3 G<H and E>F: 

a) If rs2> E — F, which requires Vg > E — F, then 

dei = [D-H- r,i]+ + [B-F- r,2]+ 

and 

mindei = [[D - H]+ + [B - F]+ - rs]+ 

b) If rs2<E-F than 

de2 = [D-H- r,i]+ + [B- r,2]+ + - r,2]+ 

and 

de2 = [B + E- 2Ti + [D- H]+ - T2]+ 

where 

Ti = min{5, E,E-F, 
72 = r/ - Ti 

For these two cases we have 

de = max{dei,de2} 



4) CASE 4 G < H and E < F. 
We have simply: 

de = [[D-H]+ + [B-F]+ -rs]+. 

For ([I4fb : 

We need to solve: 

df = min{7ii + 721 + 712 + 722} 

subj. to max{A - 711 - [F - 712]^, - 712 - [F - 712]^}^ + [B - 722 - [F - 712 
+ max{D - 722 - [H - 721]+, G - 721 - [i? - 721]^}^ < rj. 
We divide the optimization into two steps. First we solve 

dfi = min{7ii + 721 + 712 + 722} 

subj. to max{yl - 711 - [F - 712]'^, - 712 - [F - 712]+}+ + [B - 722 - [F - 712]+] 

and 

df2 = min{7ii + 721 + 712 + 722} 

subj. to max{Z) - 722 - [H - 721]^, G - 721 - [H - 721]^}^ < r/2. 

Then we solve 

df = min{d/i + d/2} 
subj. to rfi + rf2 = Vf. 

The optimization problem has the following four forms: 
1) CASE 1: E > F and G > H: 

a) rfi>E-F and r/2 > G - H, which requires rf>[E- F]+ - [G - H]+ . 
If r'f^ < max{A',B') - mm{A',B'), then 



df = dfi + df2 = [max{A', B'} - r'f^] + [D - H - r/2]" 



and 



with 



mind/i = [ma.x{A' , B') + [D - - r/]+ 
A' = [A- F -[E - F]+]^, B'=[B-F]- 



If r^.^ > max(^', B') - xma{A' , B'), then 

\A' + B' -r'fA+ 
mind/2 = 7y ^ + [D-H- 



where 



rf2 = [D-H-Rf]+ 
rfi = rf- r/2 



r'n = hi -[B- 



and 



df = maxjdji, df2} 

Jf rfi > E — F and rf2 < G — H, which requires rf > E — F. 
If < max(^', B') - min(A', B'), then 

df = [max{A', B'} - R'f^] + [D - r/2]+ + [G 

and 

mind/1 = [max(A', 5') + + G - 2ri - 

where 

Ti = min(D, G,G-H, r,) 
T2 = rs- Ti 
If > max(yl', S') - min(A', B'), then 

mind/2 = 7. ^ + [D - ii" - 



where 



and 



r/2 = min{m£ix{D, G}, r/}"*" 



r'fi = [Tfi -[B- 



df = max{d/i, d/2} 
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c) If Vfi < E — F and r/2 > G — H, which requires rj > G — H, then 

df = [E- r/i]+ + [A- r/i]+ + [B - F]+ + [D - H - r/2]++ 

thus 

df = [E + A-2Ti + [D- H]+ -T2 + [B- 

where 

Ti = min{£;, A,E - F,[rf -[B - F]+]+} 
T2 = rf-[B- F]+ - Ti 

d) If rfi<E-F and r/2 < G - il, which requires r/ < [£; - F]+ + [G - H]+, then 

d/ = [£; - r/i]+ + [A- r/i]+ + [5 - F]+ + [L> - r/2]+ + [G - r/2]+ 

Thus 

mind/ = [£; + A - 2ri + D + G - 2T2 + [B - F]+ - T3]+ 

where 

Ti = min{£;, ^, - F, [r/ - [5 - F]+]+} 

r2 = min{D, G,G-H, [r/ - Ti - [B - F]+]+} 

T3 = [r/ - Ti - T2]+ 

and hence 

df = min{(i/i, d/2, d/3, d/4} 

2) CASE 2: E > F and G < iJ: 

a) rfi > E — F which requires rfi > E — F. 

df, = [[[A - + [B- F]+ - r/i]+ + [D-H- r/2]+ 

and 

mind/ = [[[A - + [S - F]+ + [D - H]+ - r/]+ 

b) r/i < E - F. 

d/2 = [£; - r/i]+ + [A- r/i]+ + - F]+ + [D - - r/2]+ 



and 

= [[E + A- 2Ti + [D - H]+]+ - Ta + [S - 
df = max{(i/i, df2} 

3) CASE 3: E <F mdG>H. 

a) If rf2> G — H, which requires rf2> G — H, then 

dfi = [[[A - + [5 - - r/i]+ + [D-H- r/2]+ 

and 

mind/ = [[A - F]+ + [B- F]+ + [D - H]+ - r/]+ 

b) If rf2<G-H then 

df2 = [[A - + [B- F]+ - r/i]+ + [D- r/2]+ + [G - r/, 

and 

mind/2 = [[D + G- 2Ti + - F]+ + [B - F]+ - T2]+ 

where 

Ti = min{D, G,G-H, r,}T2 = r, - Ti 

and hence 

d/ = max{d/i, d/2} 

4) CASE 4: E < F andG < H. 
We simply have 

df = [[A - F]+ + [B- F]+ - r/i]+ + [D - H - r/2]+ 

and thus 

mind/ = [[yl - F]+ + [B- F]+ + [£> - - r/]+. 
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